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Probing O VI Emission in the Halos of Edge-on Spiral Galaxies 
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ABSTRACT 

We have used the Far Ultraviolet Spectroscopic Explorer to search for 
O VI AA1031.926, 1037.617 emission in the halos of the edge-on spiral galax- 
ies NGC4631 and NGC891. In NGC4631, we detected O VI in emission toward 
a soft X-ray bubble above a region containing numerous Ha arcs and filaments. 
The line-of-sight component of the motion of the O VI gas appears to match the 
underlying disk rotation. The observed O VI luminosities can account for 0.2—2% 
of the total energy input from supernovae (assuming a full O VI emitting halo) 
and yield mass flux cooling rates between 0.48 and 2.8M0yr~^ depending on 
the model used in the derivations. Based on these findings, we believe it is likely 
that we are seeing cooling, galactic fountain gas. No emission was detected from 
the halo of NGC891, a galaxy in a direction with considerably high foreground 
Galactic extinction. 

Subject headings: ISM: general — ISM: individual (NGC4631, NGC891) — 
galaxies: general — galaxies: halos — galaxies: individual (NGC 4631, NGC 891) 
— galaxies: ISM 
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1. INTRODUCTION 

In a seminal work, Spitzer (1956) predicted that the Milky Way is surrounded by a 
tenuous, hot corona. Provided that the temperature of such gas was below the escape 
temperature for a galaxy, he showed that the gas should extend several kiloparsecs above 
the midplane. Spitzer's original arguments have been extended to a more general theory 
of the dynamical evolution of the gaseous disks of spiral galaxies. This so called "galactic 
fountain" model (Shapiro & Field 1976; Bregman 1980) postulates that the hot gas produced 
by multiple overlapping supernovae would be buoyant in the thin, cold disks of gas in spiral 
galaxies. As the hot gas rises from the thin disk into the halo, it cools, perhaps forming 
condensations of neutral material which then fall back towards the thin disk as high- or 
intermediate-velocity clouds. 

The existence of such gas in the Milky Way is now well documented through observa- 
tions of strong O VI absorption (e.g. Savage et al. 2000) and X-ray emission (Burrows & 
Mendenhall 1991). The 1032/1038 A doublet transition of Li-like oxygen (O VI) is an impor- 
tant tracer of hot, collisionally ionized material. At temperatures of ~ 3 x 10^ K, where gas 
of near-solar metallicity cools most efficiently, emission from the O VI doublet can be the 
primary coolant. Hence, observations of O VI emission provide a direct probe of the cooling 
rate of material near this temperature. Furthermore, because the cooling of gas at tempera- 
tures of a few times 10^ K is so much more efficient than that of material at higher (or lower) 
temperatures, O VI emission traces rapidly cooling material that may condense into clouds 
in the galactic fountain scenario. Observations of O VI from this "transition temperature" 
gas in the halos of galaxies provide fundamental information on the circulation of material 
within spiral galaxies. 

An O VI absorption line survey of the Milky Way (Savage et al. 2003) showed that the 
observed O VI was best explained by a patchy disk corotating with the Galactic plane with 
signs of outflows. O VI observed in emission in the Milky Way has been associated with gas 
in or around the Local Bubble (Welsh et al. 2002 and references therein) and an outflow 
from the Perseus arm (Otte, Dixon, & Sankrit 2003). With our improving understanding 
of the O VI distribution inside the Milky Way, it becomes necessary to search for O VI 
emission in other galaxies as well in order to be able to distinguish between peculiarities 
and common properties of galaxies. We analyzed spectra of two galaxies (NGC4631 and 
NGC891) taken by the Far Ultraviolet Spectroscopic Explorer (FUSE). Both galaxies are 
nearly edge-on (i.e. allow clear distinction between disk and halo gas) and known for strong 
extraplanar emission in radio continuum, optical lines, and soft X-rays (SXR). Thus, they 
appear to be good candidates to search for evidence of galactic fountains. While the spectra 
of NGC891 did not reveal any O VI emission, we detected the O VI doublet in NGC4631. 
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This is the first detection of O VI in emission in a spiral galaxy other than the Milky Way 
and its companions. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. NGC 4631 

We obtained spectra with the low resolution aperture (30" x 30") of FUSE at two 
positions above the disk of NGC 4631 (see Figure 1) as part of Project P134. Position A is 
at the peak of the 1/4 keV X-ray emission detected by Wang et al. (1995) using ROSAT. 
Position B is between two Ha arcs that were believed to form the walls of an open chimney 
into the halo (Rand, Kulkarni, & Hester 1992), although a recent Hubble Space Telescope 
Ha image (Wang et al. 2001) indicates that this feature may be a superposition of several 
loops. All FUSE observations were conducted in TTAG mode where the detectors record 
the arrival time, x- and y-position, and the pulse height of each event (Sahnow et al. 2000). 
About 2/3 of each observation was conducted during orbital night. 

The O VI lines in the NGC 4631 spectra were extracted from the raw data files. Only 
events with a pulse height between 4 and 15, inclusive, were used in the analysis. Nearby O I 
airglow lines were used to determine the height of the extraction windows, and the spectra 
were summed perpendicular to the dispersion direction. Emission free regions adjacent to 
the O VI lines in the dispersion direction were used to determine the background. 

Figure 2 shows a plot of the extracted spectra for both pointings toward NGC 4631. 
The wavelength scale was determined by running the spectra through the FUSE data 
pipeline V2.0.5. Its uncertainty is about one pixel (0.007 A). Both the redshifted 1032 A 
and 1038 A line were detected on the LiFlA channel. However, the latter was blended with 
the O lAA 1039.23,1040.94,1041.69 airglow triplet. The background at wavelengths > 1040 A 
was increased due to a stripe of scattered light within the satellite. The O VI emission lines 
were also seen on the LiF2B segment for Position B; however, the lower effective area signif- 
icantly reduced the quality of the data. Therefore, we have used only data from the 1032 A 
emission line of segment LiFlA in the analysis. 

A comparison of day-plus-night versus night-only data yielded the same flux in the 
O VI A1032 line and therefore ruled out the possibihty that the emission arose from scattered 
solar O VI emission in the LiF channels. Deep integrations towards other hues of sight ruled 
out detector artifacts, scattered light, and stray light contamination at the positions of the 
detected hues. The background was not significantly higher during the sunht portion of 
the orbit; therefore, we have used the combined day and night data in our analysis. The 
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resulting signal-to-noise ratios for the integrated O VIA1032 intensities are 4.5 (Position A) 
and 6.0 (Position B). 

2.2. NGC891 

The data for NGC 891 were obtained as part of FUSE Guest Investigator Cycle 2 pro- 
gram B114. Two positions below the disk and one position above the bulge of NGC 891 
(see Figure 3) were observed. The data reduction was identical to the data reduction for 
NGC 4631 above. Approximately 85% of the observations were conducted during the night 
portion of the orbit. The observational data for both galaxies are listed in Table 1. Given 
the large random errors in the intensities of NGC 891, we do not believe that any of these 
detections are significant. 

3. RESULTS 
3.1. NGC 4631 

NGC 4631 is an Sc/SBd type galaxy at a distance of about 7.6 Mpc with an inclination 
of i 85°. The interaction between NGC 4631 and its two neighboring galaxies NGC 4656 
and NGC 4627 might have triggered the active star formation observed in NGC 4631 today. 
The galaxy possesses a complex system of Ho; filaments and bubbles (Wang et al. 2001) and 
an extended radio halo with spurs up to lOkpc above the disk (Hummel &: Dettmar 1990). 
X-ray observations of NGC 4631 revealed a SXR bubble on the north side of the disk (Wang 
et al. 1995) and a halo of gas extending up to 8kpc above the plane with temperatures of 
(2 — 7) X 10^ K (Wang et al. 2001). Both the radio and the X-ray halo appear asymmetric 
with the larger extent being on the north side of the disk. 

The two FUSE pointings for NGC 4631 were chosen to lie within the SXR bubble ob- 
served with ROSAT (Wang et al. 1995, see Figure 1). Assuming an effective area of 27.0 cm^ 
for the LiFlA segment and 900arcs^ for the aperture area, we derive O VIA1032 intensi- 
ties of 4600 ± 1000 photons s^^ cm^^ sr^^ (Position A) and 8000 ± 1300 photons s~^ cm~^ sr~^ 
(Position B). Table 2 lists the intensities (also converted to ergss^^ cm^^ sr~^) and the cor- 
responding height above the disk. The O VI A1032 intensity at Position B is a factor of five 
smaller than the 2 a upper limit of Ferguson et al. (1995), who observed the southern portion 
of the galaxy with the Hopkins Ultraviolet Telescope on the Astro-2 mission. The full extent 
of the O VI emitting gas is Tinknown. However, it is useful to consider two cases. In the first 
scenario, the O VI emission is confined to the region of the SXR bubble on the north side of 
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the disk. On the other hand, the O VI emission may extend across the full hot halo revealed 
by the Chandra observations of Wang et al. (2001). The former would correspond to a sin- 
gle large bubble breaking out into the halo, while the latter implies a full galactic fountain. 
Therefore, we will calculate the luminosity for both cases assuming that the average O VI sur- 
face brightness is equivalent to the average of Positions A and B. The O VI A1032/O VI A1038 
line ratio approaches two in optically thin gas and is reduced in optically thick gas due to 
self- absorption. We assumed the upper hmit for this line ratio of 2 in order to account for 
the O VIA1038 emission in our calculations. Ignoring extinction (see below) we calculated a 
lower hmit for the total O VI luminosity of Lovi = (6.2 ± 0.8) x 10^^ ergs s~^ for the bubble 
and Lqyi — (2.3 ± 0.3) x 10^^ ergs s~^ for the entire halo. 



3.1.1. Correction for Extinction 

The Milky Way extinction toward NGC4631 is E{B-V) = 0.017 (Schlegel, Finkbeiner, 
& Davis 1998). Using a Galactic extinction curve by Cardelli, Clayton, & Mathis (1989), we 
derived an optical depth r = 0.25, i.e. the UV flux is attenuated by a factor of 1.28. The 
corrected values for the luminosities are Lqyi = (7.9 ± 1.1) x 10^^ ergs s"^ for the bubble 
and Lqyi = (3.0 ± 0.4) x 10^^ ergs s^^ for the entire halo. The amount of extinction within 
NGC4631 is uncertain. Long slit optical spectroscopy of the Balmer lines by Martin & Kern 
(2001) has revealed evidence for significant and patchy extinction due to dust in the halo of 
NGC4631 up to 5kpc from the plane. However, another group found almost no extinction 
at a height of > 1 kpc on the northern side of the disk (Otte et al. 2001). We obtained long 
slit spectra of NGG4631 with the 2.3 m Steward Observatory on Kitt Peak, Arizona, on 
2002 April 7-9. The slit was about 4'/5 wide and orientated north-south to cross both FUSE 
positions. The measured Hq;/H/3 line ratios are comparable to those measured by Otte et 
al. (2001). We therefore assume that extinction within the halo of NGC4631 is neghgible. 



3.1.2. Comparison with X-ray Luminosity 

As mentioned earlier, the O VI emission is an important coolant for hot ionized gas. In 
fact, for gas in heating-cooling equilibrium, the O VI emission cooling at intermediate tem- 
peratures (few xlO^K) should dominate the radiative coohng in the SXR emission of 10^ K 
gas (Sutherland & Dopita 1993). The SXR intensities in the 0.1 — 2kcV range measured by 
Chandra and integrated over the area of the FUSE positions are 1.5 x 10^^ ergs s~^ cm^^ sr~^ 
(Position A) and 3.1 x 10^^ ergs s~^ cm~^ sr^^ (Position B) yielding SXR/O VI ratios of 
1.2 ±0.3 and 1.3 ±0.2 for Positions A and B, i.e. the intensities are more or less comparable 
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between the O VI doublet and the SXR emission. However, the used O VI AA1032,1038 emis- 
sion was a lower limit (see paragraph 3.1) and not corrected for Galactic extinction (because 
we did not know the Galactic extinction for the SXR emission); the SXR/0 VI ratios did 
not account for any dust in NGC4631 either. Thus, the quoted ratios are probably lower 
limits suggesting that the gas is not in equilibrium. 

The total O VI luminosity (corrected for Galactic extinction, see paragraph 3.1.1), if 
assumed to cover the whole halo of NGC4631, is by a factor of 9 lower than the two com- 
ponent fit of the X-ray luminosity (2.7 x 10^° ergs s^^, Wang et al. 1995, corrected for 
different distances). The corresponding temperature of this X-ray luminosity was approxi- 
mately 6 X 10^ K, i.e. a temperature about twice as high as the one for the peak of the O VI 
coohng curve. This again imphes that the gas is not cooling efficiently in O VI, i.e. it is not 
in equilibrium. However, the uncertainties here are higher than for the comparison above, as 
we now have compared observational data with a model using different simplifications (e.g. 
for the shape and extent of the emitting halo) . 

Wang et al. (1995) estimated that the mechanical energy input from supernovae is about 
(1.5 — 15) X 10^^ ergs s~^, i.e. the observed O VI luminosity can account for at most ~ 2% 
of the total energy input from supernovae. If the gas is in equilibrium, then the remaining 
energy input has to be used to transport the gas from the disk into the halo against the 
gravitational potential of the galaxy and against the pressure of the magnetic field lines 
and is partially lost due to cooling of ~ 10^ K hot gas by species other than oxygen. From 
the cooling flow models of Edgar & Chevalier (1986), we calculate that the full halo case is 
equivalent to a mass flux of 2.8M0yr~^ for cooling at constant density and 1.8M0yr^^ for 
cooling at constant pressure. In case of a single cooling bubble the corresponding mass fluxes 
are 0.74 Mgyr"^ (isochoric) and 0.48 Mgyr"-^ (isobaric). A rough estimate for the mass flux 
derived from i?(95'>lT' observations is ~ IMqJt:'^ (Wang et al. 1995). Despite the large 
uncertainties due to unknown extinction and the differences in the cooling models, the data 
suggest that the O VI emitting region is larger than the SXR bubble, but smaller (within a 
factor of ~ 2) than the X-ray halo observed with Chandra. This means that outflows are more 
signiflcant in the central region of the galaxy, as indicated by the fllamentary structure in the 
Ha image (Wang et al. 2001). If, in fact, there was strong, patchy extinction along our lines of 
sight, the actual O VI luminosities could be much higher than the observed luminosities and 
could account for a significant fraction of the input energy from supernovae. On the other 
hand, if the O VI emission originated in the interface between cold clouds moving through 
the hot coronal gas, the relationship between SXR and O VI emission would become more 
complicated. 
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3.1.3. Rotation Velocity 

We fitted the observed O VI A1032 emission line by a convolution of a tophat (106kms~^ 
wide to represent the filled FUSE 30" x 30" square aperture) and a Gaussian profile (to 
represent the spectral resolution of FUSE). We used calculations to determine the best 
central wavelengths and intrinsic velocity dispersions for the emission lines at Positions A 
and B. In order to estimate 1 a uncertainties for these parameters, we varied each parameter 
independently until the corresponding minimum increased by 1 relative to the best fit x^. 
The derived values are Ac = 1034.30 ± 0.09 A and FWHM = 200 ± 50kms~^ for Position A 
and Ac = 1034.37 ± 0.08 A and FWHM = 170 ±50 km s'^ for Position B. The corresponding 
radial velocities (geocentric) are v = 690 ± 30kms~^ (Position A) and 710 ± 20kms~^ 
(Position B). An offset of —9.1 and — 9.5kms~^, respectively, yields local standard of rest 
velocities. If the O VI emitting gas does not fill the LWRS aperture and is not centered 
in it, the maximum additional offset in the line centroid is 50kms~^ (for a point source 
at the border of the aperture). The line broadening due to the sateUite's orbital motion is 
< 8kms~^, i.e. well within the uncertainties of the FWHM. 

The intrinsic FWHMs, if assumed to be thermal widths, correspond to temperatures of 
(1.4 ± 0.7) X 10^ K (Position A) and (1.0 ± 0.6) x 10^ K (Position B). These temperatures 
are an order of magnitude higher than the hottest component of the SXR emission. This 
suggests that the line widths are caused by turbulent motion or differential rotation rather 
than thermal motion. From position- velocity diagrams (Rand 1994), one can expect an 
offset of about 50-100 km s~^ between the systemic velocity (606kms~^) and the velocity of 
gas corotating in the disk at the distances of the FUSE positions projected onto the major 
axis. A combination of the size of the low resolution aperture (30" or about 1.1 kpc) and 
differential rotation yields a velocity dispersion of up to 200kms~^ in the galactic plane 
near the center. The O VI velocity centroids and FWHMs closely match those of the H I 
gas in the disk of NGC4631. Thus, the hne-of-sight component of the motion of the O VI 
emitting gas appears to follow the rotation of the underlying disk gas. Due to the galaxy's 
high inclination, the velocity of the O VI emitting gas perpendicular to the disk is unknown. 
Although not conclusive, this strong relationship to the disk gas implies that the O VI 
emitting material originated in the disk. 

3.1.4. Scale Height 

If the gas has an exponential distribution in the vertical direction above the disk, we 
can calculate the density scale height of the gas from the two observations at Positions A 
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and B using: 

ln(7|) 

where and are the vertical heights of the projected aperture above the plane of the 
galaxy, and /a and Jb are the measured specific intensities at Positions A and B, respectively. 
The 2 in Equation 1 results from the fact that the intensity is proportional to the square of 
the density. Using the intensities and z heights from Table 2, we find 

/i = 8.3±4.1kpc. (2) 

The large uncertainty in the density scale height results from the large uncertainties in 
the measured specific intensities at Positions A and B. However, it appears that the scale 
height of O VI gas in NGC4631 is significantly larger than the 2.3 kpc scale height found 
by Savage et al. (2003) for the Milky Way. Using the intensities from paragraph 3.1.2, we 
find a scale height of about 6.6 kpc for the SXR emission in NGC4631, whereas Hummel, 
Beck, & Dahlem (1991) reported a scale height of about 8 kpc for the magnetic field. Scale 
heights derived from optical emission lines range between 0.9 and 1.2 kpc, not considering the 
roughly constant intensity at 2; > 3 kpc (Martin & Kern 2001). Given the uncertainties, the 
scale heights of the O VI and SXR emission and the magnetic field appear to be comparable. 



3.2. NGC891 

The Sb spiral galaxy NGG891 is at a distance of about 9.6 Mpc with an inchnation of 

i > 88.6° (Rupen et al. 1987). We observed three positions in the halo of NGC891 and did 
not detect O VI in emission to a level greater than 3a (Fig. 4). The 3cr intensity limits can 
be obtained from Table 2. The Milky Way extinction in this direction is E{B — V) = 0.065 
(Schlegel et al. 1998). Assuming a standard interstellar extinction curve (Cardelli et al. 
1989), the UV flux is attenuated by a factor of 2.4 (r = 0.88). If we assume that the O VI 
emitting halo of NGG891 is the same size as the X-ray emitting halo seen by Bregman 
& Pildis (1994), and that the average specific intensity of the halo in NGG891 is equal 
to the average at Positions 1 and 2, then the upper hmit to the O VIA1032 luminosity is 
i^i032 < 7.5 X lO-'^^ergss-^ Bregman & Houck (1997) found that the SXR (0.1-2.5 keV) 
band luminosity is 2.45 x 10^^ ergs s^^ with a cooling rate of 0.08 M0yr~^ (corrected for 
d = 9.6 Mpc). This predicts an O VIA1032 luminosity (Edgar & Chevalier 1986) of 5.7 x 
]^q37 gj^gg g-i (isochoric) and 8.8 x 10^''ergss^^ (isobaric), well within our limits. We have not 
attempted any correction for extinction in the halo of NGC 891 itself. 
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4. CONCLUSIONS 

NGC4631 has a stronger radio halo than NGC 891, but a comparable far-IR luminosity. 
Rand, Kulkarni, & Hester (1992) concluded from these observations that NGC 4631 probably 
has a larger star formation rate than NGC 891. The Hubble Space Telescope images of Wang 
et al. (2001) clearly show loops, arcs, and filaments as would be expected from galactic 
chimneys. ROSAT observations (Wang et al. 1995) revealed the presence of a large region 
of SXR emission directly above the arcs and filaments, indicating that some of the hot 
supcrnovae gas has vented into the halo. In this paper, we reported the detection of O VI 
emission at two locations above the arcs at the peaks of the SXR emission with a line-of-sight 
velocity similar to the disk gas and a FWHM broader than what might be expected from 
thermal motions of hot gas. The observation of O VI emitting gas suggests that the hot gas 
originated in the disk and is cooling. Due to the large uncertainties in the O VI intensities 
and luminosities and the range in the luminosities and mass flux rates predicted by cooling 
flow models, it is not clear whether the cooling flow in NGC 4631 is (yet) in equihbrium 
with the mass flux derived from X-ray observations. However, given all this evidence, we 
believe that, because of its edge-on nature, NGC 4631 is the flrst spiral galaxy in which O VI 
has been detected in emission unambiguously from the galactic halo and in which all the 
expected components for a full cycle of gas in a galactic fountain have been observed. The 
lack of reliable O VI measurements for NGC 891 is likely explained with the large foreground 
extinction of the Milky Way and the amount of dust in NGC 891 itself. The O VI luminosities 
predicted from cooling flow models are well within our 3(j limits for NGC 891. 

The NGC 4631 data were obtained for the Guaranteed Time Team by the NASA-CNES- 
CSA FUSE mission operated by the Johns Hopkins University. Financial support to U. S. 
participants has been provided by NASA contract NAS5-32985. The NGC 891 data were 
obtained as part of the FUSE Cycle 2 Guest Investigator Program B114. Financial support 
to U. S. participants has been provided by NASA contract NAG5-10250. JCH recognizes 
support from NASA grant NAG 5- 10957. This work made use of the NASA Extragalactic 
Database (NED) and the NASA Astrophysics Data System (ADS). 
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Table 1. OBSERVATIONS 



Position 


Program 


Right Ascension 


Declination 


Exposure Time 


Events in LiFlA 




ID 


(J2000 h m s) 


(J2000 ° ' ") 


(s) 


1032 A line 


NGC4631-A 


P1340101 


12 42 08.8 


32 34 36.0 


20954 


55± 12 


NGC4631-B 


P1340201 


12 42 08.8 


32 33 36.0 


16220 


74± 12 


NGC891-1 


B1140101 


02 22 29.0 


42 21 12.0 


31603 


15± 12 


NGC891-2 


B1140201 


02 22 40.0 


42 22 36.0 


31513 


25± 13 


NGC891-3 


B1140301 


02 22 44.8 


42 22 12.0 


18531 


12± 10 



Table 2. SPECIFIC O VIA1032 INTENSITIES 



Position 


z 


Ao32 






(kpc) 


(photons s"-*^ cm~^sr~^) (10~' 


^ergss ■'^cm ^sr ^) 


NGC4631-A 


4.8 


4600 ± 1000 


9±2 


NGC4631-B 


2.5 


8000 ± 1300 


15 ±3 


NGC891-P 


2.4 


< 2000 


< 3.8 


NGC891-2^ 


1.4 


< 2200 


< 4.2 


NGC891-3'^ 


4.1 


< 2800 


< 5.5 


^The listed 


intensity 


is a 3cr upper limit. 
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mOHl ASCENSION {-13000) 

Fig. 1. — X-ray contours superimposed on an Ha image of NGC4631 (Wang et al. 1995). 
Both FUSE positions are marked in the top panel, which shows the SXR bubble observed 
in the ROSAT 0.15-0.3 keV band. The bottom panel shows the hard X-ray contours of the 
0.5-2.0 keV band overlaid over the same Ha image as in the top panel. 
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Wavelength (A) 



Fig. 2. — FUSE day-plus-night spectra of NGC4631 at Positions A (top) and B (bottom). 
The spectra were binned by 16 pixels for display only. The O VI AA1032,1038 doublet is 
marked with arrows; airglow lines are marked with the Earth symbol. The dashed lines show 
the extraction window over which the O V1A1032 photons were counted. The continuum 
fit for this emission line is plotted as a horizontal line. Its extend outside the dashed lines 
shows the region of the spectrum used to determine the continuum. 




Fig. 3.— BV image of NGC891 (Howk & Savage 2000). The three FUSE positions are 
marked. 
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Fig. 4. — Same as Fig. 2, but for NGC891 at Positions 1-3 (from top to bottom). The 
positions of the arrows and extraction windows are based on the rotation curve by Sancisi 
& Allen (1979). 



